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ABSTRACT 

We present visible spectroscopic and albedo data of the 2.3 Gyr old Themis family and the 
<10 Myr old Beagle sub-family. The slope and albedo variations between these two families 
indicate C-complex asteroids become redder and darker in response to space weathering. Our 
observations of Themis family members confirm previously observed trends where phyllosili- 
cate absorption features are less common among small diameter objects. Similar trends in the 
albedos of large (>15 km) and small (<15 km) Themis members suggest these phyllosilicate 
feature and albedo trends result from regolith variations as a function of diameter. Observa¬ 
tions of the Beagle asteroids show a small, but notable fraction of members with phyllosilicate 
features. The presence of phyllosilicates and the dynamical association of the main-belt comet 
133P/Elst-Pizarro with the Beagle family imply the Beagle parent body was a heterogenous 
mixture of ice and aqueously altered minerals. 

Subject headings: main belt asteroids, carbonaceous asteroids, space weathering 


1. Introduction 


Space weathering studies on asteroids have pri marily focused on th e alteration of the silicate rich S- 
complex asteroids. These moderate albedo (pv '^0.22: [Mainzer et^boi 1 1 asteroids are known to spectrally 
darken, redden and have i ncreasingly supp ressed absorption bands as a function of time ( Belton et al.lll992 . 
19941 : IBinzel et al.lll99q: IChapmanlll996ll . These spectral changes are attributed to the vapor deposition 
of sub-mic ron metallic iron (SMFe) particle s onto grains during micrometeorite impacts an d solar wind 
irradiation ( YamadaetaL 1999 : Haoke 2001 : Sasaki et al. 2001 : Brunetto & Strazziilla 2005). The intrin¬ 


Mainzer et al. 


201 ih nature of carbonaceous material and the lack of prominent 


sically dark (pv ~0.06; 

absorption features at visible and n ear-IR wavelengths implied spac e weathering trends would be difficult 
to identify for C-complex asteroids (lHapkdl200ll : lMoroz et al.lll996h . However, two recent studies indicate 
significant spectral s lope variation among this class of asteroids as a function of age (iNesvorny et al.ll2005l: 
Tazzarin et al.ll200(Th . 


Principal component analyses of asteroid colors in the Sloan Digital Sky Survey (SDSS) show the 
mean slopes of C-complex asteroid families experience a decrease in slopes and become spectrally bluer 
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with age at visible wavelengths (iNesvornv et alJl2005h . In contrast, a s tudy using visible sp ectroscopic 
data from the Small Main-Belt Asteroid Spectroscopic Survey (SMASSIL IBus & Binzelll2002bh. shows the 


spect ral slopes of the C -complex asteroid pop ulation as a whole increase (redden) with age (ILazzarin et al. 


20061) ■ However, when lT,azzarin et al.l (120061) limit their analyses to the C - complex asteroid farn ilies they 
are able to reproduce the slope trends obtained by iNesvomv et alJ (120051) . iT.azzarin et al.l (120061) suggests 
the discrepancies between the two slope trends arise from a sampling effect where the average compositions 
of C-complex asteroids are not fully represented in the sampled C-complex families. 


The INesvornv et all (l2005l) and iLazzarin et al.l (|2006l) studies indicate significant evidence of space 
weathering trends among the C-complex asteroids, howev er the effects of c ompo sitional variation on these 
trends is still unclear. In the case of S-compIex asteroids, INesvornv et all (l2005l) use the Koronis asteroid 
family and Karin cluster, which formed from the breakup of a Koronis family member, to study space 
weathering tr ends while avoiding th e influence of compositional variati on. The spectral tre nds of the young 
Karin (5 Myr; Nesvomv et al.ll2002 ) and old Koronis asteroids (2.5 Gyr; Marzari et all 1995 ) confirm that the 
reddening observed among the S-complex asteroids is indeed a result of space weathering and not a product 
of mineralogical variation. The recent discovery of a sub-family o f the Themis ast eroids allows a similar 
test to be conducte d on C-complex asteroids. The Themis family (lHiravamalll918l) has ~4,000 members 


dMilani et al.ll2014) resu lting from the catastrophic break up of a 390-450 km parent body asteroid 2.3 Gyr 


ago (iMarzari et al.lll995h . A recent (< 10 Mya) break up of a Themis family memb er (D'^20 km to 65 km 
in size) resulted in the formation of the Beagle family, which contains ~60 asteroids ( Nesvornv et aklboOS ). 


These two families are the first C-complex families identified which originate from the same parent body and 
provide a unique tool for assessing C-complex space weathering trends while alleviating the mineralogical 
variations among C-complex asteroids. 

In this study we use spectroscopic and albedo data to search for space weathering trends among the 
Themis and Beagle aster oids. We compare th ese trends to those of the Veritas asteroid family, which is 
another young (8.3 Mvr: INesvornv et al.ll2003l) C-complex family. The Veritas asteroids are not related to 
the Themis and Beagle asteroids and allow us to assess whether differences in mineralogy affect the spectral 
trends in space weathering studies of C-complex asteroids. 


2. Observations 


In this paper, we present data for 52 main-belt asteroids belonging to the Themis, Beagle and Veritas 
asteroid families. Observations took place on 3.5 nights during Mar. 03, Oct. 29, Oct. 30, 2013 and Feb. 
21, 2014 at the 8.2-meter Subaru telescope on Maunakea, Hawai’i (see Tabled for log of observations). 
Spectroscopic data covering the 0.47 < A < 0.91 /im spectral range were taken with the Faint Object Camera 
and Spectrograph (FOCAS). FOCAS has two 2Kx4K CCDs with a total of 8 readout channels, each with 
512x4176 pixels dKashikawa et al.ll2002l) . The spectra were recorded in channel 3 of chip 2, which has a 
read noise of 3.4 e- and gain of 2.082 e-/ADU. We used the lowest resolution grating (75 gr/mm) and 2x2 
binning to obtain the highest signal to noise ratio (SNR) for our targets which are dominated by relatively 
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small (D <15 km) and faint asteroids. The resulting low resolution dispersion of 11.8 A/pixel is ideal 
for detecting the shallow and broad absorption features typically seen in C-complex asteroid spectra. The 
SY47 order sorting filter was used to prevent higher order contamination. Observations were taken through 
a 1” wide slit, oriented to the parallactic angle, while tracking at non-sidereal rates. Due to the lack of 
non-sidereal guiding and to prevent drifting across the slit, integration times were limited to 600 seconds. 


Family members were selected based on the dynamically derived memberships found in iNesvomv 


(l2012h . Previous stu dies show trends where the frequency of a steroids with absorption features varies as a 
function of diameter ( Florczak et al.lll999 : Fomasier et ahlboid ), so Themis and Veritas targets were limited 
to members with similar diameters as Beagle targets (D < 15 km) to avoid size related variations. Nearby 
solar analog stars were observed at airmasses similar to the asteroids throughout the night and ThAr lamp 
spectra were taken at the beginning or end of each observing run for wavelength calibration. 


3. Data Reduction and Analysis 


Data were reduced usin g the Image Reduction and Analysis Facility (IRAF) V2.14 noao longslit and 
Subaru f ocasred packages (lTodvlll986h . The reduction procedure included overscan and bias subtraction, 
image trimming, flattening, and cosmic ray removal for long exposure images. The apall package was 
used to perform background subtraction, extraction of the 


one-dimensional and associated sigma spectra from the two-dimensional images. The sigma spectrum is 
produced by measuring the sigma at each wavelength in the two-dimensional image. Following extraction, 
data were wavelength calibrated using emission lines from the ThAr lamp spectra. 


Median combined asteroid spectra were divided by nearby G2V solar analog spectra to produce re¬ 
flectance spectra for each asteroid. The solar analog stars were observed close in time to the observations 
of each asteroid and were chosen to have an airmass difference of <0.1 to enable extinction correction. The 
reliability of each solar analog star was tested by producing multiple reflectance spectra of an asteroid from 
multiple stars at similar airmasses. Only standard stars that produced consistent spectral shapes and slopes 
were used in producing reflectance spectra. Table [T] reports which solar analog stars were used to produce 
the reflectance spectra for each asteroid. Residual features due to incomplete removal of sky lines were 
removed using a median filter to search for erroneous spikes in the final reflectance spectra. 


We measured slopes using the equation for t he normalized reflectivity gradient S' (reported in units of 
%/1000 A), as defined by .lewitt & Meech (1986): 


S'iXuX2) 


fdS/dX\ 

V *^0.55 / 


( 1 ) 


where dS/dA is the slope of the reflectivity measured within the wavelength region between Ai and A 2 . 
Reflectance spectra were normalized to 0.55 /^m and spectral slopes were measured using a weighted linear 
least-squares fit to data between 0.49 < A < 0 .91 Ttm. The fits were weighted using the IRAF generated 
sigma spectrum for each asteroid. Data shorter than 0.49 /xm were not included in the fit due to potential 
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contamination from the UV absoiption below ^0-5 um created by an Fe^+ intervalence charge transfer 
atttibuted to phyllosilicates (IVilaslll994l: iHiroi et alJll996h . Asteroid slopes and 1-sigma uncertainties are 
shown in Table |2] 



phase coefficients from other C-complex asteroids does little to improve slope measurements for observa¬ 
tions at small phase angles (< 40). In addition, we do not have phase reddening coefficients for most of our 
targets, but our observations span a small range of phase angles (<20°), thus no phase reddening correction 
was applied to our reflectance spectra. 


One of the fe w features found in C -type spectra is a shallow (< 5%) 0.7 /rm absorption band created 
by phyllosilicates dVilas & Gaffevlll989l) . To search for this feature we perform continuum removal and fit 
the 0.55-0.85 /xm feature region with a 2nd or 3rd order polynomial. The continuum removed spectra were 
created by dividing the full spectrum by a linear fit to the 0.7 /im continuum shoulders (0.52-0.54 /rm and 
0.86-0.88 /xm). The continuum removed spectra are shown in Figures [Dill The band depth and center were 
determined using the reflectance value at the wavelength corresponding to the minimum of the polynomial 
fit. The band depth uncertainty was derived using the uncertainty in the sigma spectrum corresponding to 
the position of the band center and only features with band depths greater than the associated uncertainties 
were flagged as a detection. Band depths for objects where features were detected and sensitivity limits 
computed from the SNR of each spectrum are reported in Table |2] 


In addition to our Subaru data, we use visible albedos derived bv lMasiero et al.l (120111) from the Wide- 
field Infrared Survey Explorer (WISE) to characterize the albedo distributions of Themis, Beagle and Veritas 
asteroids. We used 2025, 565, and 16 Themis, Veritas and Beagle members, respectively, in our analyses 
but divided the data set into two size populations (cutoff =15 km) for consistency with our Subaru data. 
The mean values for slope and albedo, the standard error of the mean (SEM) and the standard deviations for 
each family and size range are reported in Table [3l We use the SEM to quantify the variation in our estimate 
of the mean, and the standard deviation to quantify the dispersion in our slope and albedo distributions. 


4. Results 


Members of the Themis and Beagle families show a wide range of blue to red-sloped spectra (Table 
la. However, ~90% (21 out of 23) Beagle asteroids have blue-sloped spectra, whereas 60% (13 out of 
22) of the Themis asteroids are blue-sloped. These data suggest the B eagle family is dominated by B-typ e 
asteroids, which are characterized by negatively sloped (blue) spectra (ITholenll 1984]: (Bus & Binzelll2002ah . 
The Themis asteroids show a range of C, B, E and G spec tral types (lTholenlll984h . consistent with the 
observations of other Themis members (lElorczak et al.lll999n . 
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Fig. 1.— Continuum removed spectra of individual Themis (T) asteroids, where labels report the asteroid 
ID. 


4.1, Spectral Slopes 

The slope distributions for the Themis and Beagle families are shown in Figure [H Given the larger 
fraction of blue sloped spectra, the weighted mean for the Beagle asteroids (—1.280 ± 0.003% per 1000 A) 
is significantly bluer than that of the Themis asteroids (—0.378 ± 0.003% per 1000 A). The difference in 
slopes between Themis and Beagle imply these asteroids redden by ~0.9% per 1000 A over a timespan of 
~2.3 billion years. To test the significance of fhe difference in fhe slope disfribufions of fhe fwo families, we 
applied fhe Sfudenf’s T and Kolmogorov-Smimov (KS) fesfs fo fhe Themis and Beagle dafa sefs. We find 
the difference in means is significant at the 2-sigma level. 

The Themis asteroid 18946 is extremely red (3.037 ± 0.005% per 1000 A), and is located near the 
outer region of the proper element distributions for the Themis family. Due to the possibility of 18946 being 
an interloper, we exclude 18946 from our analyses, and compute a new mean slope of —0.505 ± 0.003% 
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Fig. 2.— Continuum removed spectra of individuals asteroids. Labels show the asteroid ID as well as family 
membership, where (T) and (B) are Themis, Beagle family members, respectively. 


per 1000 A for the Themis asteroids. The reduced mean equates to a slope reddening rate of ~0.8% per 
1000 A in ~2.3 billion years. After re-applying the statistical tests, the difference in means between the 
Themis and Beagle populations is still significant at the 2-sigma level although there is a slight drop in 
significance. 

The Beagle asteroids 140399 and 8161 were each observed on two nights and have signihcantly varied 
spectral slopes (Table O. Both spectra of 140399 show the asteroid is moderately blue, however spectra of 
8161 are both blue and red sloped. Both spectra were taken during photometric conditions and produced 
using the same solar analog star, thus the large slope variation in asteroid 8161 may imply color differences 
across the surface of this asteroid. However, follow-up observations are needed to assess whether these 
trends can be reproduced. 
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Fig. 3.— Continuum removed spectra of individuals asteroids. Labels show the asteroid ID as well as family 
membership, where (B) and (V) are Beagle and Veritas family members, respectively. 


4.2. Hydration Features 


Florczak et al.l d 19990 and iFornasier et al.l(l2014t) suggest that the 0.7 /im feature is less frequent among 


asteroids with decreasi ng diameters. Data fro m the SMASS II survey show a high fraction of hydr ated 
members in the young dNesvorny et al.l 120031 8.3 Myr;) Veritas asteroid family (IBus & Binzelll2002bh . so 
we took spectra of 5 randomly selected small (D < 10 km) and 2 moderately sized (D~20 km) Veritas 
members to test if the 0.7 /um feature was observable on small diameter objects. We found 4 out of 7 (57%) 
asteroids less than 25 km in diameter show the 0.7 /um feature. Although asteroid 169282 does not pass 
our detection criteria (the depth of the band is less than the SNR limits), visual inspection of the spectrum 
(Fig. m does suggest the feature is present. Including asteroid 169282, 3 out of 5 members (~60%) with D 
< 15 km show the 0.7 /um feature, suggesting we should expect to detect the feature on small Themis and 
Beagle asteroids if phyllosilicates were present. 
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Fig. 4.— Continuum removed spectra of individual Veritas (V) asteroids, where labels report the asteroid 
ID. 


However, only one Beagle (140399) and no Themis asteroids show the feature. It is important to note 
that asteroid 140399 was imaged on two separate observing runs and only one of the spectra (Fig. ^ shows 
the absorption. Despite the lack of 0.7 /rm features in the Themis members, a small number of o bjects 
(4 of 22, or 18%) show evidence of the Fe^+ UV absorption band also attributed to phyllosilicates dVilas 


19941) . A similarly small fraction (3 of 23, or 13%) of the Beagle asteroids also show this UV feature. Due 
to the cutoff in the spectral coverage below 0.47 /im, we did not attempt to fit the Fe^+ UV feature, but its 
presence was assessed through visual inspection. Table [2] shows which asteroids were flagged as showing 
the UV feature. 
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Fig. 5.— Spectral slope gradients (as defined bv ijewitt & Meechl (Il986h ) for Themis and Beagle members. 
The weighted mean and standard error of the mean are reported. 


4.3. Albedo Variations 


The WISE survey sensitivity is strongly dependent on the distance and diameters of surveyed asteroids. 
As Beagle and Themis members have nearly identical orbital elements, the detection sensitivity in our 
analyses is limited to asteroid diameters. Therefore, due to the small size range of the Beagle asteroids, we 
limit the comparison between the three families to objects with diameters less than 15 km. Figure [6] shows 
the WISE albedo distributions as function of diameter for both the small (D < 15 km) and large (D > 15 km) 
asteroid populations in each family. The population of small Themis asteroids have a lower average albedo 
(Pv = 0.068±0.001) than the Beagle asteroids (pv = 0.079±0.005). Unexpectedly, the small diameter Veritas 
asteroids, which are similar in age to the Beagle family, have albedos (p,, = 0.069±0.001) similar to the 
older Themis asteroids within the same size range. In addition, the larger Themis asteroids (D > 15 km) are 
significantly brighter (pv = 0.075±0.001) than the smaller members. The Veritas asteroids also show some 
evidence in albedo variations between small and large members, however the difference is not statistically 
significant. 


For comparison between our results and IMasiero et al.l (1201311 . we fit gaussian profiles to the log distri¬ 
butions of the three families. The distributions in log space are shown in Figure |7] We find that our gaussian 
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Fig. 6 .— WISE albedo distributions for Beagle, Themis and Veritas asteroid families from iMasiero et al 
(l2013b . 


means are consistent with the values derived bv IMasiero et al.l (12013h . When comparing the gaussian medi¬ 
ans with the observed sample distribution, we find the gaussian values are slightly larger than the means of 
the observed population. The difference in means imply the observed population does not follow a gaussian 
distribution. The apparent non-gaussian distribution may originate from an observational b ias introduced 


by the dependence of WISE derived albedos upon previously published absolute magnitudes (IMasiero et al. 


20131). However, since the Themis, Beagle and Veritas families are located at similar heliocentric distances. 


the observational biases are similar for each family, allowing comparison of the observed albedo distribu¬ 
tions for these families. 


5. Discussion 


Spectral characteristics such as slopes and absorption features are one of the few tools available for 
assessi ng asteroid c ompositions remotely. The abunda nce of the absorber, s ize distribution of the regolith 


grain s (IClarklll999h . and space weathering processes dConel & Nashlll970l: iGold et al.lll970l : iHapke et al. 


1970ll all greatly influence fhe fealures of asferoid specfra. Allhough Ihese processes and characlerislics are 
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Fig. 7.— Log plot of WISE albedo (pv) distributions for Themis, Beagles and Veritas asteroids. Each 
histogram has been normalized to unit area for comparison. Gaussian mean and widths are reported for each 
family. 


difficult to disentangle when analyzing spectroscopic data, the sensitivity of spectroscopy to these processes 
allows it to be a powerful remote sensing tool. 


5.1. C-Complex Space Weathering Trends 


As seen in Eigure [5l the older Themis asteroids have distinctly redder slopes than the younger Beagle 
asteroids. This difference implies the spectral s lopes of C-complex a steroids become redder over time, 
which is consistent with the trends observed by lEazzarin et al.l (l2006h. Although the slope evolution of 


C-com plex asteroids appears con sistent with si l icate-r ich asteroids (INesvornv et all 120051 : lEazzarin et al. 


20061) . the albedo trends are not. iMasiero et al.l (|201 ll) find fhe younger Karin family has a lower albedo 
fhan the older Koronis family. In contrast, our data suggest the younger Beagle asteroids have a higher 
albedo than the older Themis asteroids. Euture work is necessary to assess whether the albedo trends in 
Karin and Koronis are representative of the broad S-complex weathering trends. 


The relatively young ages of the Veritas and Beagle asteroids suggest they have minimally evolved 
regoliths, yet the Veritas asteroids show a much higher fraction (57%) of hydrated members than the Beagle 
asteroids (13%). Therefore, the low fraction of features in the Beagle asteroids likely reflects an intrinsically 
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low phyllosilicate content. However, considering the higher fraction (35-60%) of large Themis members 
with phyllosilicate features (IFlorczak et al.lll999ll . we postulate the diagnostic features have been removed 
due to regolith evolution and space weathering processes. We describe this hypothesis further in § 15.21 In 
addition to the difference in phyllosilicate features, the Veritas family albedos show that they are darker 
than the Beagle asteroids, but comparable to the older Themis asteroids. Thus, these data imply a significant 
difference in the comp osition between the Veritas and Themis/Beagle families and that the trend observed by 
Nesvorny et al.l (120050 is likely a product of mineralogical variation rather than space weathering processes. 


Classical space weathering models attribute the reduced band depths , darkening and reddening of at¬ 
mosphereless bodies to the optical effects of SMFe particles (lHapkell200lh . Although C-c omplex asteroid s 
appear distinct in their mineralogy when compared to lunar soils and S-complex asteroids ( Bus et all2002 1. 
many of th e visible waveleng th absorption features in C-complex asteroids are attributed to oxidized Fe 
transitions (IRivkin et al.ll2002l . and references within). Thus, the availability of Fe and our observed spec¬ 
tral trends suggest space weathering processes on C-complex asteroids results in the production of SMFe 
particles. 

Using the age difference between the Themis and Beagle families, we estimate a slo pe increase of 
~0. 8% per 1000 A or ~0.08 /im ^ in 2.3 Gyr for C-complex asteroids. In comparison, lUazzarin et al.l 
(120060 estimates a weathering rate of 8.8 ± 4.6 x 10^^ /xm^' AU^ for the C-complex asteroid population. 


which in 2.3 Gyr equates to ^0.02iU m ^ for objects with c ircular orbits located near the Themis asteroids 


(~3.1AU). For S-complex asteroids, lUazzarin et al.l (l2006h estimates a rate of 24.9 ± 4.6 x 10 ^ /rm 


AU^, corresponding t o ^0.06 um ^ in 2.3 G yr for objects at 3.1 AU. Although no estimate was derived for 
C-complex asteroids, iNesvorny et al.l (1200511 derive a weathering rate for S-complex asteroid s of 0.01 


X log in? {t given in Myr), which equates to ~0.03 /xm ^ in 2.3 billion years. In contrast. IVernazza et al. 


(120091) suggests space weathering occurs very rapidly within the main asteroid belt, finding that S-complex 
asteroid slopes redden by as much as ~0.4 /rm^' in several million years for newly formed families. A 
direct comparison between each of these studies is difficult due to the inherent difference in each data set, 
likely causing t he significant d i fferen ce b etween each of the esti mates. However, our derived rates support 
the estimates of lUazzarin et al.l (120061) and iNesvornv et al.l (|2005|) . which suggest that space weathering is a 
long term process for both S and C-complex main belt asteroids. 


5.2. Size Related Spectral Trends 


Several groups note trends of a decreased frequency in phyllosilicate abso^tion features with de creas- 
ing size of asteroids ( Jones et al. 1990 : Vilas & Svkes 1996 : lHowell et ah 2001 : Fornasier et al. 2014 ). The 
same trend has been noted in the Themis asteroids in particular, with ~60% of astero ids greater than 50 km 
presenting phyllosilicate features in contrast to 35% for objects less than 50 km. dFlorczak et al.lll999t) . 
Additionally, our data show that the fraction reduces to 18% for Themis asteroids with diameters less than 
15 km. Furthermore, these size related trends also extend to Themis family albedos, where objects with less 
than 15 km are darker than larger asteroids. 
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Due to their lower escape velocities and the inability to retain fine particles duri ng collisions, small 
asteroids are th ought to develop relatively thin regoliths dominated by coarse grains (iMatson et alJll977l: 
Chapmanll 1978 1. Indeed, spacecraft observations of small asteroids such as Eros and Itokaw a reveal surfaces 


dominated by gravel and me ter sized boulders with fines consfricfed fo low pofenfial zones (iChapman ef al. 


20021: iMivamofo ef al.ll2007tl. These rocky su rfaces are hypofhesized fo be more resisfanf fo gardening pro¬ 
cesses dChapmanl l2002l : IWillman et al.l l2010ll . thus space weathering effects (decreased albedo and band 
depths) are likely to be more prevalent on small objects with time. Conversely, impacts upon large asteroids 
with thicker regolit hs and finer petioles may be more effective in counferacfing space weafhering fhrough 
regolifh gardening (IWillman et al.ll2010tl . thereby preserving spectral features for longer. Furthermore, lab¬ 
oratory data show that i ncreasing grain size results in a decreased reflectanc e or albedo and a reduction in 
absorption band depths ( Fischer et al.lll994 : Clarklll999 : Cloutis et ah 2011a b). We propose the size related 
absorption feature and albedo trends in the Themis asteroids originate from the combined effects of coarse 
grained regoliths and optical maturation times, which result in a pronounced reduction of absorption features 
and decreased albedos for smaller asteroids. 


In contrast to the spectral trends seen in the Themis asteroids, iBinzel et al.l (Il996t) found tha t the small 
near-E arth asteroids (NEAs) exhibit deeper absorption bands when compared to large NEAs. IChapman 
(120041) suggests these trends originate from the rapid erosion and resetting of regol iths on smaller aste roids. 
However, impact velocities are lower in the main belt than the near Earth region (iBottke et al.ll2002h . thus 
impact processes will different between the two populations and may explain the difference in spectral trends 
with size. Further observations are necessary to assess whether the albedo and absorption feature trends seen 
in the Themis asteroids are found among other main belt asteroid families. 


5.3. Aqueous Alteration in the Themis Parent Body 


Se veral main-belt comets (MBC) (iHsieh & Jewittil2006t) have been dynamically linked to the Themis 
family ( Hsieh 2009 : Novakovic et al.i2012 ). The source of sublimation in these MBCs is attributed to wate r 
ice buried beneath a layer of regolith several meters thick (ISchorghofenl2008l:IPrialnik & Rosenberg|l2009h . 
Thermal evolution models suggest the Themis/Beagle parent body was likely differentiated, with aqueous 
alteratio n taking place predominantly in the core, thus allowing the preservation of an ic y shell near the 
surface ( Castillo-Rogez & Schmidt 2010 : Grimm & McSween 1989 : Cohen & Coker 2000). 


Due to the young age of the Beagle family and the small fraction of hydrated Beagle members, we 
believe that the low frequency of features represents an intrinsically low phyllosilicate content in the Beagle 
parent body. The detection of phyllosilicates in the Beagle asteroids has important implications for geo¬ 
physical models of the Themi s/Beagle parent body, considering the MBC 133P/Elst-Pizarro is a possible 
member of the Beagle family (INesvomv et al.ll2008h . If indeed a member, then the presence of I33P and 
its volatiles implies the Beagle parent body originated from an ice-rich layer in Themis/Beagle parent body. 
However, the small but notable fraction of members with the UV and 0.7 /rm features suggest the Beagle 
parent body was not homogeneous, but rather a mixture of ice and aqueously altered material. 
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We propose two possibilities for heterogeneity of the Beagle parent body. The first is that heterogene¬ 
ity ma y re present non-uniform a lteration of the Themis/Beagle parent body. Conversely, iRivkin & Emery 
(1201011 and ICamoins et al.l (l2010ll have potentia l ly det ected water ice features on the largest member of the 
family, 24 Themis. ICastillo-Rosez & Schmidn (l2010ll suggest 24 Themis is the remnant silicate rich core 
onto which components of its icy crust may have reaccreted following the collisional disruption of the par¬ 
ent body. It is possible the Beagle parent body may be an object which reaccreted non-native components, 
resulting in the present heterogeneity of the Beagle asteroids. Exploring the viability of these two scenarios 
is beyond the scope of this paper, however coupled geophysical and dynamical modeling of the two families 
may provide valuable insights into the evolution of this body. 


6. Conclusions 

Through comparison of visible spectra and albedo data of the Themis, Beagle and Veritas asteroids, we 
find the following: 

• Space weathering results in an increase of spectral slopes and a decrease in albedo among C-complex 
asteroids. 

• The trends in decreased phyllosilicate features and albedos with decreasing diameter likely result from 
variations in regolith properties as a function of age and diameter. 

• The notable fraction of Beagle members with phyllosilicate features indicate the Beagle parent was a 
heterogeneous mixture of ice and aqueously alterated minerals. 

• The apparent mineralogical differences between the Veritas family and the Themis and Beagle fami¬ 
lies highlight the importance of accounting for mineralogy when interpreting space weathering trends 
across the broad population of C-complex asteroids. 


We would like to say mahalo nui loa to Jan Kleyna, Takashi Hattori and Bin Yang for their help and 
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Table 1. Observing Geometry and Conditions 


Asteroid 

Vmag 

a 

e 

i 

rt 

A* 
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No. of 

Airmass 

Solar 

Airmass 

ID 


(AU) 


(deg) 

(AU) 

(AU) 

(deg) 

Run 

Spectra 

(Object) 

Analog 

(Stand) 
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17.5 
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0.157 
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3 

2 
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4 
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Table 1—Continued 
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4 

2 

1.33 

HD95868 

1.28 

Veritas 

5592 

17.3 

3.173 

0.062 

8.511 

3.254 

2.696 

15.777 

3 

2 

1.05 

HD73708 

1.03 

7231 

17.1 

3.168 

0.075 

9.429 

2.947 

2.215 

14.862 

3 

2 

1.20 

HD220764 

1.28 

8624 

18.9 

3.166 

0.032 



2.384 

11.364 

4 

2 

1.04 

HD42160 

1.14 

62512 

19.0 

3.165 

0.065 




3.871 

4 

2 

1.05 

HD73708 

1.03 

82961 

18.6 

3.163 

0.057 

8.317 

3.145 

2.163 

3.138 

3 

2 

1.09 

HD19061 

1.02 

169282 

20.0 

3.166 

0.082 



2.536 


4 

2 

1.04 

HD257880 

1.02 

215270 

19.5 

3.172 

0.079 


2.951 



2 

2 

1.15 

HD19061 

1.20 


Note. — Orbital elements (a, e, i) ephemerides, and observational data for each asteroid. Observing run numbers correspond to Run 1= Mar. 
02, 2013, Run 2 = Oct. 29, 2013, Run 3 = Oct. 30, 2013, Run 4= Feb. 21, 2014. '^heliocentric distance ^geocentric distance ^phase angle 
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Table 2. Subaru Asteroid Characteristics 


Asteroid 

ID 

Hmag 

(km) 

Slope 

(%/ioooA) 

SNR* 

UV^ 

Abs 

0.7 /xm 
Abs 

0.7 /xm Band 
Detection Limit (%) 

Obs Run'' 

No. 


Themis 

7511 

13.7 

7.09 

-0.896 ± 0.020 

158.5 

Y 

N 

0.6 

3 

8697 

13.5 

9.16 

-2.165 ±0.008 

25.0 

N 

N 

4.0 

4 

15572 

13.4 

8.12 

0.088 ± 0.004 

83.6 

Y 

N 

1.2 

3 

15772 

13.4 

8.82 

-1.917 ±0.007 

40.4 

N 

N 

2.5 

4 

16877 

13.6 

9.19 

-0.065 ± 0.010 

62.5 

N 

N 

1.6 

4 

18946 

14.4 

6.36 

3.037 ± 0.005 

37.0 

N 

N 

2.7 

1 

26633 

14.1 

8.06 

-0.025 ± 0.004 

95.9 

Y 

N 

1.0 

3 

28789 

14.4 

6.36 

-0.900 ± 0.006 

66.0 

N 

N 

1.5 

3 

31093 

13.4 

9.43 

0.293 ± 0.013 

82.6 

N 

N 

1.2 

4 

31334 

13.6 

10.18 

0.631 ±0.019 

101.5 

N 

N 

1.0 

4 

38365 

14.1 

7.79 

2.209 ± 0.017 

35.4 

N 

N 

2.8 

1 

47920 

13.4 

11.25 

0.437 ± 0.012 

129.7 

Y 

N 

0.8 

3 

51143 

13.4 

9.93 

-1.083 ±0.007 

79.0 

N 

N 

1.3 

3 

51579 

14.2 

6.01 

0.807 ± 0.004 

33.1 

N 

N 

3.0 

3 

54320 

14.4 

im 

-1.072 ±0.015 

124.0 

N 

N 

0.8 

3 

58231 

14.0 

7.64 

-1.426 ±0.012 

126.5 

N 

N 

0.8 

3 

58302 

14.2 

6.91 

-1.473 ±0.004 

71.9 

N 

N 

1.4 

1 

82646 

14.2 

6.97 

0.115 ±0.007 

50.5 

N 

N 

2.0 

4 

91965 

14.1 

6.47 

0.250 ± 0.005 

55.2 

N 

N 

1.8 

1 

123646 

14.6 

5.80 

-2.557 ± 0.005 

45.8 

N 

N 

2.2 

1 

132845 

15.3 

4.24 

-0.688 ± 0.007 

21.2 

N 

N 

4.7 

3 

149084 

15.0 

4.78 

-1.161 ±0.004 

27.1 

N 

N 

3.7 

1 

Beagle 

7968 

15.6 

4.87 

-1.132 ±0.004 

13.5 

N 

N 

7.4 

3 

8161 

13.3 

10.51 

-1.787 ±0.004 

99.3 

N 

N 

1.0 

2 




0.964 ± 0.010 

31.2 

N 

N 

3.2 

4 

19923 

13.8 

8.08 

-0.128 ±0.008 

85.3 

N 

N 

1.2 

4 

21569 

13.9 

13.65 

0.455 ± 0.015 

55.5 

N 

N 

1.8 

1 

26140 

14.1 

7.83 

0.418 ± 0.004 

54.0 

N 

N 

1.9 

3 

55188 

13.6 

8.69 

-1.924 ±0.002 

32.9 

N 

N 

3.0 

2 

63951 

15.0 

5.54 

-1.664 ±0.003 

39.6 

N 

N 

2.5 

3 

71489 

13.8 

8.13 

-0.242 ± 0.003 

28.2 

N 

N 

3.5 

3 

84230 

15.1 

5.49 

-2.572 ± 0.005 

32.2 

N 

N 

3.1 

3 

88079 

14.9 

5.49 

-3.353 ± 0.005 

35.6 

N 

N 

2.8 

1 

104671 

15.4 

4.36 

-1.261 ±0.003 

55.0 

N 

N 

1.8 

2 

116557 

15.2 

4.30 

-0.810 ± 0.005 

33.0 

N 

N 

3.0 

3 

140399 

14.6 

6.39 

-1.555 ±0.003 

40.7 

N 

N 

2.5 

3 




-2.340 ± 0.012 

55.8 

N 

3.4 ± 1.3 

1.8 

4 

143633 

15.7 

3.53 

-0.423 ±0.011 

14.9 

N 

N 

6.7 

4 
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Table 2—Continued 


Asteroid 

Hmag 

D" 

Slope 

SNR* 

UW 

0.7 /xm 

0.7 jim Band 

Obs Run*^ 

ID 


(km) 

(%/ioooA) 


Abs 

Abs 

Detection Limit (%) 

No. 

144732 

15.6 

3.70 

-1.082 ±0.003 

17.9 

N 

N 

5.6 

2 

150625 

16.1 

3.26 

-2.444 ± 0.009 

22.0 

N 

N 

4.5 

1 

166780 

15.7 

4.09 

-2.506 ± 0.004 

58.5 

N 

N 

1.7 

2 

178208 

15.6 

4.34 

-1.895 ±0.008 

41.6 

N 

N 

2.4 

4 

180762 

15.5 

3.29 

-1.379 ±0.013 

10.8 

N 

N 

9.3 

1 

182228 

15.6 

4.34 

-0.710 ± 0.008 

23.4 

N 

N 

4.3 

1 

202360 

15.6 

3.57 

-0.469 ± 0.003 

21.3 

N 

N 

4.7 

2 

203061 

16.3 

3.02 

-0.363 ± 0.007 

26.3 

N 

N 

3.8 

3 

211009 

15.7 

3.98 

-2.058 ± 0.010 

70.9 

N 

N 

1.4 

4 

Veritas 

5592 

11.7 

22.09 

-1.366 ±0.011 

111.4 

N 

3.0 ± 0.7 

0.9 

3 

7231 

12.2 

18.33 

-0.588 ± 0.007 

103.9 

Y 

1.1 ±0.8 

1.0 

3 

8624 

13.6 

9.50 

0.167 ±0.002 

49.9 

N 

1.9 ± 1.8 

1.9 

4 

62512 

14.7 

5.76 

1.211 ±0.013 

73.8 

Y 

1.4 ± 1.1 

1.3 

4 

82961 

14.1 

6.93 

-2.582 ± 0.010 

65.7 

N 

3.0 ± 1.2 

1.5 

3 

169282 

14.6 

6.27 

-0.704 ± 0.006 

26.9 

N 

N 

3.7 

4 

215270 

15.1 

5.34 

1.411 ±0.002 

30.2 

N 

N 

3.3 

2 


Note. — Spectral data of asteroids observed with Subaru telescope. ° Diameters (D) were computed using mean family 
albedos when diameter measurements were not available from WISE. * Mean signal to noise ratios, measured using data 
between 0.49 < A < 0.91 jum. Indication where absorption is seen in the spectrum either in the UV or at 0.7 /im. If 
visible, the 0.7 flm band depth is reported. Detection limits derived from the mean signal to noise ratio of each spectrum. 
® Observing runs 1, 2, 3, and 4 correspond Mar. 02, 2013, Oct. 29, 2013, Oct. 30, 2013, and Feb. 21, 2014 respectively. 


Table 3. Asteroid Family Characteristics 


Family 

Slope 

%/ioooi 

^slope 

Albedo 
(< 15 km) 

^albedo 

Albedo 
(> 15 km) 

^albedo 

Beagle 

-1.280 ± 0.003 

1.080 

0.0794 ± 0.0045 

0.0179 

_ 

— 

Themis 

-0.378 ± 0.003 

1.342 

0.0680 ± 0.0006 

0.0236 

0.0747 ± 0.0013 

0.0214 

Veritas 

0.196 ±0.001 

1.532 

0.0695 ± 0.0010 

0.0245 

0.0619 ± 0.0064 

0.0255 


Note. — Slopes re present weighted me ans for each family, and albedo means were derived from 
reported WISE values i Masiero et al.ll20ll ). The uncertainty reported with the means are the standard 


error of the mean. The sigma values represent the standard deviation and show the dispersion for each 
population. 













